Hubbard model with orbital degeneracy and integer or noninteger filling 
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The photo-emission spectrum and the specific heat are studied in the Hubbard model of small 
clusters including orbital degeneracy. We focus on the degeneracy and valence dependence in the 
limit of a large Coulomb interaction. For integer valence, it is found that the degeneracy increases 
the width of the photoemission spectrum and it reduces the contribution from the charge degrees of 
freedom to the specific heat. A deviation from integer valence reduces the width of the photoemission 
spectrum. 
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I. INTRODUCTION 

Ni metal provides a typical example of the pho- 
toemission spectrum from a relatively strongly corre- 
lated system.tl The width of the Ni 3d bandcl is re- 
duced by about 25 % compared with one-particle (den- 
sity functional) calculations .0 As a further signature of 
many-bodis effects, the spectrum has a satellite ("6 eV 
satellite" )p which is displaced by roughly the Coulomb, 
interaction U between two 3d electrons on a Ni atom.El 
In Ni the many-body effects are primarily due to the 
localized nature of the 3c? orbital. The 3d band has a 
noninteger occupation (~ 0.6 holes), since the 4s orbital 
is also partly occupied. 

A different example of a rather strongly correlated sys- 
tem is provided by alkali-doped Ceo compounds, A3C60 
(A= K, Rb). It is believed, that the alkali atoms are 
almost completely ionized,ErEI and the three alkali elec- 
trons are therefore donated into the t\ u orbital of the Ceo 
molecule, which has a threefold orbital degeneracy. Since 
in this case essentially only one orbital (t\ u ) is partly oc- 
cupied, it has an almost integer occupancy. These sys- 
tems have slight deviations from stoichiometry, but the 
occupancy should nevertheless be almost integer. This 
is in contrast to the noninteger occupancy of the local 
orbital in many other correlated systems. 

For K3C60 it has not been possible to determine the 
dispersion of the bands using angular resolved photoe- 
mission, due to the presence of orientational disorder and 
problems related to the angular resolution in relation to 
the small size of the Brillouin zone. Measurements of the 
specific heat show, however, that there is no enhancement 
of the specific heat beyond what-iKpuld be expected from 
the electron-phonon interactionEt] This is quite surpris- 
ing in view of the expected strong correlation effects in 
these systems. 

The on-site Coulomb energy U hetween two electrons 
on the same Ceo molecule is largefl compared with the 
one-electran_width W of the t% u band, with U /W ~ 
1.5 — 2.5.EJO For such a large ratio one might have ex-. 
pected the system to be a Mott-Hubbard insulator.EJ 
We have shown, however, that at half-filling the ratio 



of U/W where the transition takes place grows with the 
orbital degejuepjey N of the t\ u orbital as %/iV or some- 
what £asterc3'c3 In the large U limit it was shown that 
an extra electron in a system with otherwise integer oc- 
cupancy, e.g., half- filled, can move more efficiently if the 
orbital degeneracy N is large. The reason is that the 
extra occupancy can move through a hop of any of .the 
electrons on the site with the extra occupancy.E2nL3 In 
the half-filled case, of the order N electrons can hop to 
a neighboring site. For a diffeceiit integer filling than 
half-filling, this effect is reducedllil 

We observe that in a strongly correlated system hop- 
ping is in general reduced compared with the correspond- 
ing noninteracting system. Due to the arguments given 
above, the degeneracy nevertheless enhances the effects 
of hopping on the Mott transition compared with what 
might have been expected based on the importance of 
hopping for the noninteracting system. It is then inter- 
esting to ask if this may also happen for other properties. 

Here we study photoemission spectroscopy (PES) for 
the Hubbard model of some finite systems. We show 
that for these systems and integer occupancy, the width 
of the PES spectrum is larger than the noninteracting 
band width W for a large U and N > 1, except for fill- 
ing one and 2N. These are therefore examples where 
correlation effects increase the band width. The width is 
reduced as U is reduced or as the valence starts to deviate 
from an integer and may then become smaller than W. 
Such systems with noninteger occupation are therefore 
cases where correlation effects reduce the band width. 
Systems with noninteger occupancy also develop a satel- 
lite which moves roughly linearly with the value of U, 
while such satellites have very small weights for systems 
with integer occupancy. We have given a brief discussion 
of the angular integrated photoemission spectrum for a 
system with integer valence before. t3 

In addition we discuss the contribution to the specific 
heat from charge degrees of freedom. We show that the 
orbital degeneracy can reduce this contribution, which 
may be the reason for the lack of an electronic enhance- 
ment of the specific heat in K^Ceo- 

After giving the formalism (II. A) and presenting some 
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general arguments (II. B), we show results for the photoc- 
mission spectrum for a diatomic molecule (Sec. II. C), a 
linear chain (Sec. II. D), a six atom cluster (Sec. II. E), a 
four molecule model of K3C60 (Sec. II. F) and with mul- 
tiplet effects (Sec. II. G). The specific heat is treated in 
Sec. Ill and the results are discussed in Sec. IV. 



II. PHOTOEMISSION SPECTRUM 



(5) 



The extra occupancy on site 1 can now hop to the neigh- 
boring sites i without any extra cost in Coulomb energy. 
We can form a state where the extra occupancy is on the 
neighboring site i 



A. Formalism 

Below we study multi-band Hubbard models 

<ij>ma 

i {ma)<(m' a') 



(1) 



where i, m and a are site, orbital and spin indices, re- 
spectively. The orbital index can take N values. The 
hopping integrals are given by t im Jm < and U is the on- 
site Coulomb integral. All multiplet effects have been 
neglected, except in Sec. II. G. The creation operator for 
an electron with the quantum numbers i, m and a is 
given by i\)\ ma - We are interested in the spectrum for 
removing an electron with a given quantum number k 
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The spectrum is then given by 

Pk{e) = 

J2\(M-l,n\^k\M,0)\ 2 S(e- 



E (M) 



(2) 



(3) 

E n (M-l)), 



where \M, 0) is the initial ground-state with M electrons 
and \M — l,n) is an excited final state. The correspond- 
ing energies are Eq(M) and E n (M — 1), respectively. We 
use Lanczos method to calculate the ground-state and 
the spectrum.E3 In most, but not all, cases we assume 
that the hopping matrix elements are "diagonal" in m 
and m 



t; 



t5 mm ', if i and j nearest neighbors; 
0, otherwise. 



(4) 



Here we assume t to be negative. 



B. General considerations 

We first address the hopping of an electron or a hole 
in a system which otherwise has integer occupancy. We 
consider the large t/-limit. We construct a Neel state |0) 
with N electrons of the same spin per site. We now add 
an extra electron to site I and obtain the state 



(0) 



We have here taken into account that any of the N spin 
up electrons on site 1 can hop to site i, as is illustrated in 
Fig. 0a. Assuming hopping matrijt,elements of the form 
(0), we obtain the matrix elemcntEJ 



(7) 



(i\H\l) = VNt, 



where t is the one-particle matrix element. Thus the 
hopping matrix element has been enhanced by a factor 
of y/N. In this argument we have assumed that the elec- 
tron hops against a background of a classical Neel state. 
For certain systems, we have constructed states where 
the hoppiag energy of the extra electron or hole is more 
negative.lla The hopping may therefore be enhanced by a 
somewhat larges-jjactor than y/N in the large U case for 
certain systems.t2l 

In Fig. |l|b we show schematically a doped system. The 
electron added to the central site in an inverse photoc- 
mission process sees a neighboring site which also has 
an extra electron. This additional electron blocks one 
hopping channel. In the large U case the effect is, how- 
ever, much more dramatic, since hopping to a site with 
an extra electron is then completely suppressed. 

The doping of the system has a strong effect on the 
satellite structure. Let us consider the large U limit for a 
system with the average number K + x electrons per site. 
Then a fraction 1 — x of the sites have K electrons and 
a fraction x have K + 1 electrons. If an electron is now 
removed from a site with K electrons, there is a strong 
coupling to final states with K — 1 electrons on that site. 
The energy of such a state is about U higher than for a 
final state where all sites have K or K + 1 electrons. The 
result,is a satellite at about the energy U below the main 
peakE The weight of this satellite is reduced as 1 — x is 
reduced. 
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P±(e) = a±6{s-E (2)±t), 



(9) 



a) 



b) 



FIG. 1. Hopping possibilities in the undoped (a) and 
doped (b) system. In the undoped system the extra occu- 
pancy can hop in N ways to neighboring sites. In the doped 
case hopping to other sites with an extra occupancy is blocked 
if the Coulomb interaction is large. 
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FIG. 2. Photoemission from a diatomic molecule with or- 
bital degeneracy N and an integer (TV) number of electrons 
per site at the "Brillouin zone centre" (+) and "boundary" 
(-). The dispersion (band width) E+ — E- and the ratio 
W-/w+ of the weights of the two peaks are shown. W is the 
splitting (band width) between the bonding and antibonding 
levels for the noninteracting system. The figure illustrates 
how the dispersion increases with the Coulomb interaction U 
and becomes proportional to N for large U. 



C. Diatomic molecule 



where a± = (1 + x 2 ± ^/l + x 2 )/[2{l + x 2 )} and E (2) = 
2\t\(x - VI + x 2 ), with x = U/(2\t\). For U -> the 
weight a- of the antibonding state goes to zero, and for 
U = the spectrum has only one peak. For U > the 
separation between the peaks for p + and p„ is 2\t\ = W, 
where W is the one-particle band width. Thus the width 
of the angular integrated, PES spectrum is the full one- 
particle band width W,V3 and not the width of the occu- 
pied part. We can think of this width as a dispcrsional 
width, since the two peaks in the angular integrated spec- 
trum correspond to two different values of k. In the solid, 
these k values correspond to the Brillouin zone centre and 
boundary, respectively. 

We next consider the two-site model for an orbital de- 
generacy N > 1. We first consider the case with 2N 
electrons, i.e., an integer number of electrons per site. 
For the model (Eq.(|])) of the hopping used here, the 
one-particle band width remains W = 2\t\ independently 
of N. Fig. shows the splitting E + — E- of the peaks in 
p±(e) as a function of U for different values of N. The 
splitting grows with U from its noninteracting value W 
to its large [/-limit NW. In the large [/-limit the weights 
w± of the peaks in p± become equal. As shown in the 
Fig. H this limit is reached quite slowly as U/W grows. 
In addition to these peaks, there are satellites at about 
2U lower energy. These satellites have, however, very 
small weights. The reason for the large splitting between 
the bonding and antibonding peak is the degeneracy N. 
In the case of a two-site model it can be showr£3 that the 
hopping is enhanced by a factor N instead of the factor 
^fN in the simple argument above. 
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We first consider the case of a diatomic molecule with 
only two sites and with two elections. This model has 
been treated by Harris and LangetO in the absence of or- 
bital degeneracy. An electron is removed from a bonding 
or an antibonding orbital for the orbital quantum number 
771 and a =|, 



(8) 



FIG. 3. Photoemission from a diatomic molecule with or- 
bital degeneracy N = 2 and a noninteger (2.5) number of 
electrons per site. The dispersion (band width) E- — E+, 
the satellite to main peak splitting Em — Esat and the ratio 
of the satellite and main peak weights wsat/wM are shown. 
W is the splitting (band width) between the bonding and 
antibonding levels for the noninteracting system. The fig- 
ure illustrates how the dispersion is reduced as the Coulomb 
interaction U is increased. 



and the result is averaged over m. For the case N = 1 
the problem can easily be solved analytically, and we find 



We next consider the two-site model with 2N + 1 elec- 
trons, i.e., a noninteger number (N + i) of electrons per 
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site. In the large [/-limit, the ground-state now has iV+1 
electrons on one site and ./V electrons on the other site. 
If an electron is removed from the site with N + 1 elec- 
trons in a photoemission process, the system is left with 
N electrons on each site. As discussed above, hopping 
is now completely suppressed in the large U limit, and 
we expect a negligible dispersion in the final state. If, on 
the other hand, the electron is removed from the site with 
N electrons, the result is a satellite, as discussed above. 
These results are illustrated in Fig. || for the case N = 2. 
For U = the electron can be removed from a bonding 
(k = +) or antibonding (k = — ) orbital, giving a disper- 
sion W. In this limit p± have just one peak each. As U 
is increased, the dispersion (separation between the main 
peaks E±) decreases and goes to zero for U — > oo. At the 
same time satellites develop at an energy about U below 
the main peaks. The weights of these satellites approach 
the weight of the main peaks in the large [/-limit. Since 
there are several satellites with slightly different energies, 
we have calculated the average energy of the satellites in 
the two spectra p± . The splitting between the main peak 
and the average satellite position was then averaged for 
the two spectra using the total satellite weight in each 
spectra as a weight factor. 

For the diatomic molecule with N > 1 we thus find 
that the dispersional band width increases with U for in- 
teger occupancy (N) but decreases with U for half inte- 
ger occupancy (N + 1/2) relatively to the noninteracting 
width W. 



D. Linear chain 



The one-dimensional Hubbard model has been studied 
very extensively in the literatureJlj We first consider a 
linear chain with six sites, periodic boundary conditions 
and N = 1. Fig. |Ja shows the angular resolved spectrum 
for different values of k. The figure illustrates that there 
is a substantial dispersion as in the diatomic molecule. 
Thus the width of the occupied part of the spectrum is 
only 0.25IF for noninteracting electrons, but for inter- 
acting electrons (U /W = 30) the width is about Wt3 for 
N = 1. The spectrum is more complicated than for the 
diatomic molecule, having several peaks for each fc-valuc. 
In the large [/-limit, appropriate here, the spectrum can 
be interpreted in terms of holons and spinonsO 

Fig. shows the angular resolved spectrum for 

N = 2. The spectrum has more structures than for 
N = 1 (Fig. ^a) but the gross features are very simi- 
lar. A major difference is, however, that the energy scale 
is about a factor of two larger. Thus the width of the 
spectrum is now almost twice (~ 1.8) the one-particle 
band width W. This increased width is not due to the 
creation of new satellites with a large binding energy, but 
due to a rather uniform expansion of the energy scale for 
all the features present for N = 1. This is illustrated by 
the centre of gravity which has a stronger fc-dependence 



for N = 2. Thus we may think of the increased band 
width as to a large extent being due to an increased dis- 
persion. In this sense the situation is rather similar to 
the case of the diatomic molecule considered before. 
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FIG. 4. Photoemission spectrum of a six atom linear chain 
with six electrons (integer occupancy=l, doping 8 — 0) as a 
function of k and e. The orbital degeneracy is N = 1 in a) 
and -/V = 2 in b). The one-particle band width is W and 
U/W = 30 is very large. The arrows show the centre of 
gravity of the main part of the spectrum, excluding the con- 
tribution from any satellite at about e/W ~ —30 or lower. 
For plotting reasons a Gaussian broadening of 0.04VF (full 
width half maximum) has been introduced. The energy zero 
is given by the leading peak. Comparison of a) and b) shows 
that the energy scale of the spectrum has been expanded by 
almost a factor of two due to the degeneracy. 
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FIG. 5. Photoemission spectrum of a six atom linear chain 
as a function of doping S, where the number of electrons is 
(2+<5)6. The orbital degeneracy is N = 2 and the one-particle 
band width is W . The figure illustrates how the width of the 
spectrum is reduced when S is changed from 4/6 to 1/6 or 
from to -4/6. Observe the dramatic difference between dop- 
ing 5 = 1/6 and S = —1/6. For integer fillings 1 (5 = —6/6) 
and 2 (5 — 0) and 3 (8 = 6/6) the spectrum is broad. 

We next consider the results for a doped system with 
(2 + 6)6 electrons and N = 2. The angular integrated 
spectra are shown in Fig. ||. As the doping 6 becomes 
increasingly negative, the spectrum becomes narrower. 
The same tendency is seen when 6 is positive and reduced 
from 4/6 to 1/6. For 5 — 1/6, the width of the spec- 
trum is essentially zero (apart from an artifical broad- 
ening used for plotting reasons). The explanation is the 
same as for the diatomic molecule. With doping 6=1/6 
the system has two electrons on five sites and three elec- 
trons on one site. To obtain a contribution to the main 
peak the electron must be removed from the site with 
three electrons. We then obtain a final state with two 
electrons on each site. Since hopping then costs the en- 
ergy U, it is strongly suppressed and the width of the 
spectrum goes to zero. This general tendency also shows 
up for intermediate dopings, and it reflects the reduced 
hopping possibilities for the photoemission hole if it is 



surrounded by other holes. It is instructive to notice the 
asymmetry between doping 6 = 1/6 and 6 = —1/6. In 
the first case hopping is completely suppressed in the 
final states corresponding to the leading peaks, as dis- 
cussed above. In the case 6 = —1/6, on the other hand, 
the probability that the photoemission hole can hop to a 
neighboring site would only be reduced from 1 to 5/6, if 
the occupancies were spatially uncorrelated, and there is 
a small reduction in the width. In inverse photoemission 
the situation is reversed, with a large reduction in the 
width for 6 = —1/6 and a small reduction for 6 = 1/6. 
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FIG. 6. The same as in Fig. ||b but for the doping 
5 = —2/6 (ten electrons). Comparison with 5 = (Fig. 
^>) shows that the spectra with the lowest centre of gravity 
(k = 2-7r/3, 7r) have lost weight and that the centre of gravity 
of each spectra tends to move upwards. 

To see how this happens in more detail, we show in Fig. 
^the angular resolved spectrum for 6 = —2/6. Compared 
with the 6 = results in Fig. ||b, the spectrum for k = ir 
has a small weight. Since the centre of gravity of this 
spectrum is low, the reduction of its weight leads to a 
reduced width of the angular integrated spectrum. The 
same applies to a smaller extent to the k = 2ir/3 spec- 
trum. Furthermore, the centres of gravity of the spectra 
for k = 7r/3 and 27r/3 have moved upwards. This can be 
viewed as a reduction of the dispersion, and it also con- 
tributes to the reduced width of the angular integrated 
spectrum. This upward movement of the centre of grav- 
ity is partly due to an upward movement of certain peaks 
and partly due to a reduced weight of peaks with large 
binding energies. 

Fig. [7] shows the angular resolved spectrum for a pos- 
itive doping. In this case the centres of gravity have 
moved up strongly for k = 0,7r/3 and 27r/3. In addi- 
tion, the spectrum for k = ir now has a negligible weight, 
while for 6 = it contributed substantially to the width 
of the spectrum due to its low centre of gravity. In this 
case the reduction of the width of the angular integrated 
spectrum is, however, mainly due to a reduced dispersion, 
while the situation is less well-defined in Fig. 0. 
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FIG. 7. The same as in Fig. [ija but for the doping 5 = 2/6 
(fourteen electrons). Comparison with 5 = (Fig. Wb) shows 
that the centres of gravity of the spectra for k = 0,7r/3 and 
27r/3 have moved upwards, strongly reducing the dispersion. 
Furthermore, the spectrum with the lowest centre of gravity 
(k — 27r/3, 7r) has lost most of its weight. 

It is also interesting to consider integer filling away 
from half-filling. The arguments in Sec. II. B cap_be 
straightforwardly generalized to any integer filling.Efl If 
the filling is K < N, there are K electrons which can 
fill a hole on the site where a hole was created in the 
photoemission process. If the filling is K > N there are 
2N — (K — 1) holes which can be filled on the site where 
a hole was created in the photoemission process. Above 
we found that the width of the photoemission spectrum 
grows linearly with the number of hopping possibilities 
for a linear chain. For the fillings K = 1 and K = 3 
this suggests the band widths W and 2W, respectively, 
for N = 2. These results are confirmed by the explicit 
calculations shown in Fig. In inverse photoemission 
we expect K + 1 hopping possibilities for K < N and 
2N — K hopping possibilities for K > N. 

Up to now the linear chain has been discussed in the 
limit of U/W » 1. Fig. § illustrates the validity of 
these conclusions for intermediate values of U/W. The 
figure shows the width 



W = 2^S 2 - SI 2 , 



(10) 



where S% is the ith moment of the spectrum. Also shown 
is the centre of gravity relative to the leading peak. The 
figure illustrates that the spectrum is broadened even for 
small values of U. For U/W > 1 the degeneracy further- 
more becomes important and the width is substantially 
larger for TV = 2 than N = 1. 

Fig. H shows the spectra for U /W — 1.5 for N = 1 (a) 
and N = 2 (b). The figure illustrates that a substantial 
part of the large width in the N = 2 case comes from 
the spectra for k — 27r/3 and k = 7T, which have small 
weights but structures at large binding energies. Such 
contributions would probably be very hard to detect ex- 
perimentally. There is, however, also a downward shift 
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FIG. 8. Width and the negative of the centre of gravity 
(relative to the leading peak) of the spectrum for N = 1 and 
TV = 2 as a function of U /W. The horizontal arrows show the 
results for U/W = 30. The figure illustrates how the effects 
discussed remain important down to U/W ~ 1. 



a) 

k=0 


Ji 


N=l 
8=0 

UAV=1.5 


k=jt/3 


£ 


i 

L 


k=2Tt/3 


1 

A 

A 1 — " — 




k=7t 


1 

. . A . 


-1 r- 



-2 -1.5 -1 -0.5 0.5 

E/W 



b) 
k=0 



N=2 
8=0 

U/W=1.5 



k=ic/3 



k=27[/3 



^ 1 a — + = h 



k=rc 



-1.5 -1 -0.5 0.5 

E/W 



FIG. 9. Same as in Fig. § but for U/W = 1.5. The figure 
illustrates that the effects in Fig. ^ are substantially reduced 
at intermediate values of U, but that some effects still remain. 
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FIG. 10. The same as in Fig. but for doping 3 = 8/6 
(0.67 holes per site) and U/W = 0.8 (a) and U/W = 1.0 (b). 
These parameters may be typical for Ni metal. Compari- 
son of (a) and (b) illustrates how the peaks above ~ —W/2 
move upwards as U is increased ("band narrowing") while 
the peaks below ~ — W move downwards ("satellites"). For 
(7 = the width of the photoemission spectrum is 0.75W. 



Finally we consider parameters which may be consid- 
ered as representative for Ni metal. Thus we have chosen 
U/W = 0.8 (Fig. 0a) and U/W = 1.0 (Fig. [[ojb) and the 
doping S — 8/6, which corresponds to 0.67 holes per site. 
According to Auger measurements, U/W for Ni is close to 
but slightly smaller than oneBj The number of 3d-holes 
in Ni is about 0.6. As discussed earlier, because of the 
noninteger number of electrons per sites, we expect satel- 
lites that move down in energy with U. To distinguish 
between these types of satellites and the "main" band, 
we have performed the calculation for two values of U. In 
Fig. |h] the peaks below about —W/2 move downwards 
as U is increased, while the peaks above —W/2 move up- 
wards. Although U/W is not large in this case, we can 
therefore nevertheless identify the peaks below — W/2 as 
satellites. This means that the main band has only a 
width of about O^SI'F. This should be compared with 
the width 0.75W of the the occupied part of the band 
for [7 = 0. This means that the band width has been 
reduced by about 40 %. This is a bit larger than the ex- 



perimental result (25%), but it illustrates that the model 
gives a qualitatively correct band width, in spite of its 
extreme simplicity. 

To understand the band narrowing, we observe that 
for a site with two electrons the probability that a given 
neighboring site also has two electrons would be 1/6 if the 
sites with two electrons were uncorrelated. There is then 
one hopping possibility for the hole to such a site, while 
hopping to all other site costs the energy U. This would 
then suggest a reduction of the band width to W/6. Ac- 
tual calculations for U/W — 30 give reduction by a factor 
somewhat smaller than five. The much less dramatic re- 
duction of the band width in Fig. |l^ is due to the value 
of U being intermediate. 
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FIG. 11. Photoemission spectrum of a cluster with six 
atoms in a rectangular arrangement (2 x 3) with six electrons 
(integer occupancy) and U /W = 30 for N = 1 a) and N — 2 
b). The arrows show the centre of gravity of the main part of 
the spectrum, excluding any satellite at e/W ~ —30 or lower. 
Comparison of a) and b) shows how the energy scale has been 
expanded by about a factor of 1.5 for N = 2 compared with 
N = 1. 
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E. Six atom "rectangular" cluster 

As an further example we study a cluster with six 
atoms arranged in a rectangle with three times two sites, 
using periodic boundary conditions. Fig. |ll| a shows the 
spectrum for the orbital degeneracy N = 1. In this case 
the width of the photoemission spectrum is evea-larger 
than the one-particle band width W. The reasonE3 is the 
frustration in the six atom cluster due to the periodic 
boundary conditions. Because of this, the one-particle 
band width W = 5|t| is smaller than one might have ex- 
pected (6|i|) from the hopping integral and the number 
of nearest neighbors (=3). The frustration effect on the 
many-body spectrum for a large U /W is smaller, giving a 
width of the spectrum of about 5.75 \t\ = 1.1 5W already 
for N — 1 . The width of the angular integrated spectrum 
is to a substantial amount due to the dispersion. 

Fig. |lT| b shows the spectrum for N = 2. The spec- 
trum shows large similarities with the case N = 1 (Fig. 
|il"la) , except that the energy scale is larger by a factor of 
about 1.5. This behavior lies in between y/~N and N and 
has been discussed in Rcf. [Hj Fig. [ll] illustrates that 
the increase of the width of the spectrum with N is not 
limited to the linear chain but a general phenomenon. 



F. A 3 C 6 o (A= K, Rb) 

The photoemission from A3C60 (A=K, Rb) is inter- 
esting, since the experimental width of the spectrum is 
about 1.5 eV,E!J while the occupied width of thc.-pne- 
particle spectrum is only of the order of 0.3 eV.Q As 
discussed in the introduction, however, the experimen- 
tal spectrum is angular resolved, and it is not known if 
this large width is due to dispersion or satellites. We have 
considered multi-band Hubbard model in Eq. (|l|) to de- 
scribe the Ceo ti u level with degeneracy three. The hop- 
ping matdx. elements are obtained from a tight-binding 
theory,E§&E3 which takes the orientational disorder in 
A3C60 into account. We use parameters which give a 
one-particle full band width W = 0.6 eV. The angular 
integrated photoemission spectrum is shown in Fig. |lj. 
The total width of the spectrum is about 1.5 eV. This 
width is larger than y/NW ~ 1 eV. This may be due to 
the width growing somewhat faster than y/~N and to the 
system being frustrated because of the fee lattice. 

The calculated width of the spectrum is of the same 
order as observed experimentally. In contrast to the ex- 
perimental spectrum there is, however, no weight in the 
range -0.3 to -0.7 eV, while the experimental spectrum 
is smoothly reduced with increasing binding energy. The 
calculated spectrum also has much too much weight at 
small binding energies. The spectrum in Fig. [L2| can 
therefore not explain experiment. The main contribu- 
tion to the broad experimental spectrum is therefore 
probably from-phonon and plasmon satellites, as sug- 
gested earlier .Ell The degeneracy mechanism discussed 



here should, however, make some contribution to the 
weight at large binding energies Hj 



CL 
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FIG. 12. Photoemission spectrum from a model of four 
Ceo molecules for U = 1.3 eV. The figure illustrates that 
apart from the main peak there is additional weight at about 
-1 to -1.5 eV. This weight does not move as U is increased 
and it is related to the three-fold orbital degeneracy. 
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FIG. 13. The dispersion C+ — CL and the width 
2-y/ S2 — Sf of the photoemission spectrum as a function of 
the multiplet splitting 4K for a two-site system, where K is 
the multiplet integral. C± is the centre of gravity for k = ±. 
A large (U /W = 30) and an intermediate (U jW = 1.5) value 
of U were considered. 



G. Multiplet effects 

An interesting question is how these results are influ- 
enced by multiplet effects. If a site with an extra electron 
or hole is in a given multiplet state, the hopping of the 
extra occupancy to a neighboring site with a similar en- 
ergy may be restricted, because the multiplet effects lift 
degeneracies. 
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Here we consider a model, where the multiplet effects 
are described by the exchange integral K between two 
different orbitals and the difference SU = U xx — U xy be- 
tween the direct Coulomb integral for equal and unequal 
orbitals. Here we use SU = 2K. Thus we add the multi- 
plet terms 



Hi -<>''^J //„„•//,,„. 2^ 2^ 



to the Hamiltonian in Eq. (0). This part of the Hamil- 
tonian has multiplets with the energies 0, 2K and 4K 
relative to the lowest multiplet, i.e., a splitting of AK . 

We consider a system with two sites, the orbital degen- 
eracy N = 2 and four electrons, i.e., in teger oc cupancy. 
Fig. |H shows results for the width 2^/ S% — S'f and the 
separation (C+ — C_) in the centres of gravity of the spec- 
trum, where Si is the ith moment. For a large U (Fig. 
|i"3|a) the width of the spectrum is gradually reduced with 
the multiplet splitting as expected, but the reduction is 
moderate. The " dispersion" measured through C+ — C_ 
drops is reduced much more. The reason for this drop is 
a change in the character of the spectrum. For K = 
the spectra for k = ± have essentially just one peak each. 
As K is increased, new satellites are formed on the en- 
ergy scale W (not K\), which lead to rapid reduction 
in the difference C+ — C_, without changing the total 
width of the spectrum drastically. The creation of new 
satellites are related to changes in the ground-state wave 
function. Due to the integer occupancy in the ground- 
state, the hopping energy is very small, and the energy 
scale over which the multiplet effects become important 
is very small. Fig. |l3| b shows results for an intermedi- 
ate value of U. In this case the hopping in the initial 
state is more important, and a small value of K has a 
small influence on the ground-state. Both C+ — C_ and 
2 \J S2 — S 2 are now reduced with AK at a much slower 
rate. Due to the smaller value of U, however, both quan- 
tities are smaller for K = than in Fig. [l3|a. The results 
above suggest that multiplet couplings reduce the effects 
discussed in this paper, but that these effects are still 
present for intermediate multiplet splittings. 



III. SPECIFIC HEAT 

If the dispersion of the states is indeed increased due to 
the degeneracy, this should also show up as a reduction 
in the specific heat. We thus consider the specific heat 



dE{T) 
dT ' 



where E(T) is the energy of the system as a function 
of the temperature T. To perform the calculation at 
finite T, we use a technique developed by Jaklic and 
Prelovsek,c3 and based on the Lanczos method. In this 
technique a certain number of states Nq are chosen ran- 
domly, and for each state \n) the appropriate expectation 
values (n\exp(—H/T)\n) and (n\exp(—H/T)H\n) are cal- 
culated using the Lanczos method. We here use the grand 
canonical ensemble, i.e., the number of electrons is not 
fixed. In this approach the specific heat is independent of 
the degeneracy for noninteracting electrons, as it should. 
The use of the canonical ensemble leads to a moderate 
degeneracy dependence even for U = 0, due to the fi- 
nite size of the system and the change of the number of 
electrons with degeneracy, i.e., the system being closer to 
infinite for a larger N. 

We consider the Hubbard model for a six-atom cluster 
with the atoms arranged in a 3 x 2 rectangle, as above. 
The hopping integrals in Eq. (^) are used. In the large 
tZ-limit the total energy at half-filling for a system with 
M sites and the orbital degeneracy N is 



1 



E(T = 0) = -MN(N - 1)U. 



(13) 



Here we have neglected small terms of the order —t 2 /U. 
For a very large temperature we obtain 



E(T 



(^) 2 MN(2N- 



1)U. 



The energy increase per electron is then 

^[E(T = 00) -E(T = 0)] = \U, 



(14) 



(15) 



i.e., independent of the degeneracy N. Integrating 
c v / (MN) over all temperatures should then give a quan- 
tity which is independent of N for large U. In this 
limit, we therefore expect the specific heat to be reduced 
with N but to remain appreciable over a correspondingly 
larger temperature range. 

We may think of the contributions to the specific heat 
as resulting from the excitation of spin and charge de- 
grees of freedom. Due to its finite size, the system has 
a band gap even for small values of U. For U/W =1.5 
and N = 1 the gap is E g = 0.76W and for N = 2 we find 
Eg = 0A5W. For T much smaller than the band gap, 
charge degrees of freedom are not excited, and only spin 
excitations contribute to the specific heat. Fig. [l4| shows 
the spin correlation function 



1 



ANM 



(ij)m 



im j ml 



(16) 



(12) 



where the summation is over nearest neighbors and af m 
is the Pauli spin matrix a z for orbital m on site i. The 
figure illustrates that the spins become uncorrelated over 
an energy scale of about a few tenths of W. This energy 
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scale is smaller than the band gap and for the corre- 
sponding temperatures the spin excitations should dom- 
inate the specific heat. The specific heat is also shown 
in Fig. [L4 . The present system may be too small to 
describe the spin degrees of freedom of a large system, 
since at low temperatures only a very small number of 
states contribute to the specific heat, and since we do 
not know if these few states are representative of a larger 
system. Thus the sharp peak in c v for small T"s may be 
a defect of the small size of the system. We are inter- 
ested in metallic systems, i.e., systems where the band 
gap goes to zero with increasing system size. For these 
systems the charge degrees of freedom should become in- 
creasingly important for low temperatures as the size of 
the system grows. Because of this and because of the 
uncertainty in the contribution from the spin degrees of 
freedom we here concentrate on the charge degrees of 
freedom, dominating for T > 0.2IT for the six atom clus- 
ter. 

For low temperatures in this range, the specific heat 
is smaller for N — 2 than for N = 1, in agreement with 
the expectations above. The specific heat decreases more 
slowly with T for N = 2, remaining appreciable up to 
larger values of T. This is required to satisfy Eq. (16). 
Finally we remark that although we have not discussed 
the spin degrees of freedom further, these degrees give 
an important contribution for many strongly correlated 
systems. 
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FIG. 14. Specific heat for a six-atom 3x2 cluster with 
U/W = 1.5 as a function of temperature T and degeneracy 
N. The result for U — is also shown. The spin correlation 
function {erf a*-) (for N = 1) illustrates that the spin degrees 
of freedom contribute only for small T, while the charge de- 
gree of freedoms contribute for larger T. The quantities have 
been normalized to the number of electrons NM. 



IV. CONCLUDING REMARKS 



We have studied the effects of orbital degeneracy N 
and doping on the photoemission spectrum in the limit of 
a large Coulomb interaction U . We find that the degener- 
acy increases the dispersional width of the spectrum for 
integer occupancy. Deviations from integer occupancy 
tends to reduce the width of the spectrum, in particular 
if the occupation is slightly larger than integer. These 
effects are gradually reduced as U is reduced or as mul- 
tiple! effects are taken into account. We also find that 
the degeneracy tends to reduce the contribution to the 
specific heat from the charge degrees of freedom, as one 
would expect from the increased dispersion. 

An interesting approach for treating, correlation effects 
is the dynamical mean field theory.E3 In this approach 
the self-energy is assumed to have no wave vector de- 
pendence, jwhich is correct for a system with infinite 
dimension. cj'EZI The present results show that in a Hub- 
bard model with integer occupancy, the orbital degener- 
acy increases the k dependence for intermediate to large 
values of the Coulomb interaction U. This suggests that 
the dimension where the dynamical mean field theory 
becomes accurate increases with the orbital degeneracy. 



It should be emphasized that we have here studied a 
(degenerate) one-band Hubbard model. For each cor- 
related system one has to consider whether or not this 
model is appropriate before applying the results above. 
For certain systems, e.g, heavy fermion and High T c com- 
pounds, the interaction between different types of or- 
bitals plays an important role, and it would be interesting 
to study how this influences the considerations above. 
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